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ABSTRACT. 

Our  results  show  that  a  generalized  model  that  includes  propagation  delays  ^d 
signal  decorrelation  can  reproduce  numerous  observed  features  in  the  spectral  modulation 
patterns  of  regional  phases  from  quarry  blasts.  Analyses  of  spectra  associated  with  various 
regional  arrivals  of  the  same  event  revealed  differences  in  modulation  stren^  and 
apparent  shifts  in  modulation  patterns  These  may  be  associated  with  source  finiteness, 
propagation  delays  and  signal  decorrelation.  Potentially  these  features  can  be  us^  to  derive 
information  with  regards  to  the  spatial  extent  and  layout  of  shot  arrays.  Modulation  patterns 
or  the  locations  of  first  minima  in  the  spectra  are  generally  not  determined  solely  by  the 
tinip.  duration  of  the  shot  sequences  as  often  claimed,  but  are  complex  functions  of 
numerous  factors.  Lack  of  spectral  modulation  does  not  necessarily  indicate  a  single 
explosion  or  the  absence  of  ripple  firing  either.  Since  this  project  has  just  begun,  the  results 
presented  here  are  relevant  only  a  to  a  subset  of  the  tasks  to  be  performed  under  this 
contract 
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OBJECTIVES. 


The  objective  of  this  project  is  the  in^lementation  of  various  advanced  time  series 
processing  techniques  for  small  array  data  and  their  implementation  in  Geotool  and, 
ultimately,  various  AFTAC  systems  for  seismic  monitoring.  The  project  includes  the 
impl^entation  of  site-specific  automatic  phase  recognition  algodtiims,  sig^  enhancement 
methods,  and  deconvolution  algorithms  combined  with  signal  complexity  discrimination. 


RESEARCH  ACCOMPLISHED. 

Spectral  modulation  (scalloping)  is  an  effective  discriminant  for  distinguishing 
ordinary  ripple-fired  quarry  blasts  fiom  large,  single  explosions  and  earthquakes.  It  cannot 
be  used,  however,  to  distinguish  between  single  chemical  explosions  and  small  nuclear 
blasts,  neither  can  it  be  used  discriminate  between  ripple-fired  quarry  blasts  and  multiple 
large  explosions,  some  of  which  could  be  nuclear,  llius  as  a  discriminant,  scalloping  has 
limitations.  Nevertheless  spectral  scalloping  can  provide  information  about  the  spatio- 
temporal  stracture  of  sources. 

Data  Analyses. 

The  regional  phase  spectra  presented  in  this  paper  were  produced  by  a  uniform 
procedure.  For  a  given  event  all  traces  were  inspect^  and  only  the  sensor  ouq>uts  with 
uncompromised  data  quality  were  kept  All  events  analyzed  also  had  very  high  signal-to- 
noise  ratios  with  the  even  weakest  arrivals  (such  as  Pn),  high  above  the  background  noise 
level.  The  various  regional  arrivals  were  windowed  using  the  same  chosen  window  length 
of  12.8  seconds  for  each  type  of  arrival.  The  Parzen  tapering  fimction  was  used  on  each 
window  which  heavily  tapered  off  both  ends  of  each  window.  The  position  of  the 
windows  were  chosen  such  that  the  arrival  to  be  analyzed  was  close  to  the  center  of  the 
window  in  order  avoid  reduction  of  its  amplitude  by  the  taper.  AftCT  Fourier  transforming 
and  squaring  the  amplitude  spectra  of  each  sensor  for  a  given  arrival  the  restilting  power 
spectra  were  smooAed  by  a  5  point  moving  average  and  averaged  over  all  the  chosen 
sensors  of  the  array. 

The  most  common  modulation  pattern  observed  in  regional  seismograms  is 
nearly  identical  for  the  various  phases  in  the  seismogram  (Figure  1).  Such  modulations 
can  be  easily  seen  in  spectra  computed  for  various  time  windows  or  continuous 
sonograms.  When  the  causative  time  patterns  of  shots  are  simple,  such  as  a  pair  of 
collocated  explosions,  we  also  can  discern  the  time  pattern  in  cepstral  analyses,  some 
seismograms  there  are  some  finer  observable  ^ctral  features  ^at  can  be  explained  by 
source  finiteness  effects.  These  observable  effects  include  phase-velocity  dependent 
frequency  shifts  (cff  rescaling)  in  the  spectral  modulation  patterns  (Fi^re  2)  and 
differences  in  the  degree  of  modulation  for  the  various  regional  seismic  arrivals  (Figure 
3). 


Theoretical  Modelling  of  Spectral  Modulation. 

A  common  assumption  made  in  various  studies  of  ^ectral  modulation  is  that  the 
modulation  patterns,  especially  the  frequency  of  the  first  spectral  miniTnum,  are  simple 
functions  of  source  duration,  i.e.  the  time  length  of  the  shot  sequence.  It  is  easy  to  show 
that  this  is  not  the  case.  The  time  delays  between  shots  at  quarries  are  generally  much 
shorter  than  the  digitizing  time  interval  at  most  sensors  of  regional  arrays  (25  msec).  Thus 
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the  arrays  essentially  "see"  a  simplified  time  function  which  can  be  derived  by  smoothing 
the  original  time  sequence  with  a  sine  function  and  resampling  it  at  a  25  msec  sampling 
rate.  Ctae  can  easily  verify  that  this  procedure  does  not  change  the  resulting  spectra  below 
the  Nyquist  frequency  of  20  Hz.  By  generating  shot  time  sequences  that  have  the  s^e 
duration  but  have  different  variations  in  charge  sizes  with  time  we  show  that  the  locations 
of  the  first  spectral  minimum  can  be  shifted  significantly.  It  is  also  easy  to  design  ripple 
firing  sequences  with  approximately  Gaussian  envelopes  when  smoothed  and  resampl^  at 
40  Hz  that  would  not  display  any  modulation  despite  tiieir  finite  duration.  Thus  the  spectral 
modulation  patterns  are  not  sin:Q)le  functions  of  d^tion. 

Obviously  the  model  of  linear  superposition  of  shots  with  identical  waveforms 
cannot  fiiUy  explain  the  observations  of  d^erences  in  modulation  of  the  various  regional 
arrivals  presented  above.  What  is  needed  is  a  model  that  incorporates  the  effects  of  si^al 
decorrelation  and  propagation  delays.  Descriptions  of  firing  practices  at  some  mines 
indicate  the  use  of  very  large  shot  arrays  with  dimensions  close  to  a  kilometer  and  firing 
time  delays  of  the  order  of  seconds  (Kim  et  al  1994).  When  a  shot  array  is  even 
moderately  large,  for  instance,  when  its  dimensions  are  of  the  order  of  hundreds  of  meters 
we  can  no  longer  assume  that  the  waveforms  of  two  widely  separated  shots  are  the  same. 
The  waveforms  originated  by  the  various  charges  are  bound  to  be  decorrelated,  by 
invoking  reciprocity  arguments,  in  a  way  similar  to  waveform  decorrelation  observed  at 
small  arrays.  The  fi^uency-limited  correlation  coefficients  (coherences)  of  waveforms  at 
small  arrays  were  found  to  be  approximately  exponentially  dependent  on  the  product  of 
distance  between  two  observation  points  and  fipequency. 

Analyses  of  waveform  decorrelation  at  NORESS  indicated  that  the  waveform 
decorrelation  is  stronger  for  for  phases  with  low  propagation  velocities  such  as  Sn  and  Lg 
compared  to  fast  travelling  phases  such  as  Pn  and  Pg.  The  stronger  decorrelation  for  Sn 
and  Lg  may  also  be  related  to  the  more  complex  mode-ray  structures  of  these  phases  (Der 
and  Baumgardt  1994).  There  are  also  indications  that  the  rate  of  decorrelation  may  depend 
on  the  direction  of  the  wave  propagation  relative  to  the  direction  chosen  in  the  horizontal 
plane  (Der  et  al  1988).  For  simplicity  in  this  paper  we  shall  ignore  this  directionality  and 

will  assume  a  decorrelation  factor  proportional  to  exp(-afd)  ,  where  a  is  a  constant 
determining  the  rate  of  decorrelation  /  is  frequency  and  d  is  the  intersensor  diatance.  The 
rate  of  waveform  decorrelation  will  obviously  depend  on  the  degree  of  heterogeneities  in 
the  local  structure  and  path  and  thus  will  vary  from  place  to  place.  The  second  factor  that 
must  be  considered  for  large  shot  arrays  is  the  propagation  delay  between  pairs  of  charges 
which  for  shot  patterns  of  the  order  of  hundred  of  meters  are  not  negligible  compared  to 
the  absolute  time  lags  among  shots.  These  delays  will  be  added  to  Ae  absolute  time 
differences  and  produce  essentially  a  Doppler  effect  For  nonsymmetrical  space-time 
distributions  of  shots  thus  will  give  a  directional  dependence  to  the  spectra  of  quarry  blasts 
observed.  Moreover,  since  the  additive  relative  propagation  delays  d^end  on  the  dominant 
phase  veocity  of  each  arrivaL  the  spectral  modulations  will  not  be  the  same  for  all  phases. 
The  formula  that  approximately  describes  the  observed  power  q)ectra  including  all  these 
factors  is 


Y(ct))  =  X^\ Xj^co ) p  +  2  ^ \Xj(a) ) ) CicD,A y^) | cos (oiTj-r^  +  d jj^ 


(1) 
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where  the  factor  C(g>^)  is  the  coherence  between  two  waveforms  corresponding  to  the 

angular  frequency  (o  and  the  distance  A,  djk  is  the  propagation  delay  for  a  phase  between 
the  two  sensors  j  and  k  which  can  be  computed  using  the  average,  dominant  phase 
velocity  of  the  phase  in  question  (such  as  Pn,  Pg,  Sn  or  Lg).  The  formula  states  that  the 
total  observed  spectrum  is  made  up  of  the  individual  spectra  of  the  sources  plus  all  the 
cross-spectral  contributions  between  all  the  possible  distinct  pairs  of  sensors.  Given  the 
fact  that  according  to  experience  spectra  of  individual  e3q)losions  are  generally  smooth, 
almost  all  of  the  modulation  is  cau^  by  the  cross  spectrd  contributions,  Le.  fre  second 
sum  on  the  tight  side  of  (1). 

In  our  simulations  we  have  assumed,  somewhat  arbitrarily,  the  dominant  phases 
velocities  of  8.2,  6.0, 4.7  and  4.0  for  Pn,  Pg,  Sn  and  Lg  respectively.  We  have  varied  the 
decorrelation  rates  to  study  their  effects  for  various  types  of  geological  situations.  Seismic 
reciprocity  as  applied  to  regional  signal  decorrelation  measurements  at  "small"  arrays  of  a 
few  km  diameter  implies  that  waveforms  observed  at  a  remote  sensor  from  similar  sources 
displaced  over  small  distances  will  be  similarly  decorrelated.  In  orda:  to  optimize  signal 
processing  at  regional  arrays  signal  decorrelation  was  measured  at  the  sites  of  small  arrays 
in  Scandinavia  ^ykkeltveit  et  aL  1983).  The  results  produced  indicated  that  the  degree  of 
decorrelation  increased  in  the  order  of  Pn-Pg-Sn-Lg  and  that  it  also  increased  with 
increasing  frequency.  In  this  paper  we  shall  assume  a  decorrelation  function  of  the 

direction-independent  form  of  exp(-o^A),  where  a  is  a  constant,/  is  frequency  and  A  is 

the  distance.  The  costant  a  will  undoubtedly  depend  on  the  local  geological  heterogeneity 
of  the  particular  area  where  the  sensor  or  shot  array  is  located  and  thus  will  be  varied. 

Applying  Equation  (1)  to  a  simulated  linear  array  demonstates  that  the  underlying 
model  can  reproduce  both  the  relative  displacement  of  the  nulls  and  maxima  of  the 
modulation  patterns  and  the  differences  in  the  strength  (peak-valley  contrast)  of  modulation 
for  the  various  phases  (Figure  4).  The  shifts  in  modulation  patterns  for  this  case  are 
lunambiguous  arid  azimuth-dependent 

A  more  realistic  areal  pattern  we  have  simulated  consisted  of  5  rows  of  ten  shots 
which  the  firing  in  the  successive  rows  overlapped  (Figure  5).  Patterns  of  this  kind  were 
described  by  Chapman  et  al  (1991).  Similarly  to  observed  data,  differences  in  the 
frequencies  of  various  features  in  the  modulation  patterns  are  small  in  this  example,  but  the 
strength  of  modulation  decreases  strongly  from  ^  to  Lg.  The  two-dimensiond  layout  of 
this  shot  pattern  makes  the  Doppler  shifo  conplex  and  hard  to  see.  This  figure  provides  an 
explanation  why  shifts  in  modulation  patterns  are  rarely  seen  in  real  observations  although 
most  shot  patterns  are  generally  large  enough  to  produce  significant  propagation  delays 
compared  to  the  absolute  time  delays  in  a  pattern. 

RECOMMENDATIONS  AND  FUTURE  PLANS. 

Obviously  spectral  modulation  patterns  in  regional  arrivals  contain  information 
about  die  spatio-ten^xaral  layout  of  the  shot  arrays  used.  This  information  can  potentially  be 
used,  m  conjunction  of  independent  information,  to  assess  the  sizes  of  shot  arrays  and  in 
discrimination.  Observation  of  modulation  patterns  from  various  directions  frrom  the  source 
may  enable  us  to  infer  tire  spatial  layout  of  the  shotarrays.  Since  this  inversion  constitutes  a 
complex  nonlinear  estimation  problem  combinatorial  optimization  methods  (such  as  genetic 
algorithms)  may  be  used  advantageously  in  such  applications. 

We  find  this  part  of  our  investigation  nearly  complete,  during  the  future  work  on 
this  project  we  shall  concentrate  on  the  other  tasks  describ^  above. 
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